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Abstract 
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fluorescence by the atoms may be completely quenched, as a result of quantum interference, and that 
when the atoms' detuning is smaller than their linewidths, extremely narrow fluorescence features emerge. 
Furthermore, the system emits a two-photon bound state which can display spatial oscillations/quantum 
beats, and can be tuned from bunched to anti-bunched statistics as the total photon energy is varied. 
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I. INTRODUCTION 



There has been substantial recent interest in the study of photon-atom interactions, where 
microwave [1-5] and optical [6-17] photons are confined to a single-mode waveguide. From a 
practical point of view, quantum states of light are important carriers of information in quantum 
information and quantum computing systems. The use of waveguides to connect qubits can enable 
entanglement transfer [16, 18], and is important for integration. From a more basic point of view, 
the one-dimensional nature of photon states in a single-mode waveguide leads to a number of 
novel physics effects in photon-atom interactions, as well as device possibilities. For example, 
confinement of photons to one dimension enables the complete reflection of a single photon from a 
two-level atom [2], and the full inversion of an atom with a single-photon pulse [10]. Furthermore, 
it has been shown that when two photons scatter off a two-level system, a photon-photon bound 
state and an associated background fluorescence emerge [3]. Recently, logic operations at the single 
photon level through the use of a three-level system in a waveguide have been investigated [19]. 

In this work, we consider a waveguide coupled to two non-identical quantum two-level systems, 
as shown in Fig. 1, and study few-photon transport. By two-level atoms, we are primarily concerned 
with on-chip atomlike objects, such as superconducting qubits [1] or quantum dots [7]. These 
objects can be described by the same two-level Hamiltonian as a real atom, hence our nomenclature. 
However, in contrast to real atoms, one of the distinct properties of these atomlike objects is their 
tunability [5, 20-22]. Our predicted effects exploit the tunability of these atomlike objects. 

Here we show that the two-atom system can be solved exactly in the two-photon Hilbert 
space, using input-output formalism [23] adapted for the calculations of few-photon Fock-state 
transport [24]. The results point to a rich set of physics, some of which may be important to 
device applications. Previously, this system was studied at the one -photon level [8, 25], and was 
shown to exhibit a single -photon transmission spectrum that is the analogous to electromagnetically 
induced transparency (FIT) [8]. Quantum entanglement and modification of spectral features via a 
resonant laser [26], as well as spatial modulation of spontaneous emission decay [27] were both 
studied in two identical two-level atoms in free space. Multiphoton scattering in a waveguide was 
studied only in multiple identical atoms [28]. Here we show that allowing for non-identical atoms 
enables fluorescence linewidth narrowing and quenching, as well as new capabilities to design and 
control the properties of photon-photon bound states. These two-photon bound states can exhibit 
bunching or anti-bunching statistics. The two photons forming the bound states can have very 
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FIG. 1 . (Color online) Schematic representation of a pair of nonidentical two-level atoms in a waveguide 
geometry. ^12 and pi 2 denote the incoming and scattered photons, respectively. Qi 2 and ti^ are the 
transition frequencies and decay times of the two atoms, respectively. 

different frequencies. Moreover, the properties of the bound state, including its spatial extent, are 
strongly dependent on the resonant frequencies of the atoms. None of these characteristics have 
been observed in waveguide-atom systems consisting of either a single two-level atom or a single 
three-level atom. 

Two-photon bound states represent a composite particle of photons, and are of substantial 
interest in quantum lithography and imaging [29]. Our results show that in the two-atom system 
there is enhanced capability for generation and control of such a composite quantum object. Also, 
one typically expects two-photon bound states to arise from effective photon-photon attraction, and 
thus one typically expects the two-photon bound state to exhibit a bunching behavior. Indeed, only 
bunched two-photon bound states have been seen in all waveguide-atom [3, 30], and waveguide- 
nonlinear-cavity [14] systems previously considered. In this context, our result, showing an anti- 
bunched two-photon state, is counter-intuitive, and points to the substantial richness in physics of 
photon-photon interaction in the two-atom system that is qualitatively different from all previously 
considered systems. 

II. SYSTEM HAMILTONIAN AND EQUATIONS OF MOTION 

For the system shown in Fig. 1, the waveguide supports both left and right propagating photon 
modes, though the photon-atom interaction is entirely contained in the even subspace [3]. This 
subspace features a chiral photonic band interacting with two atoms, which in the rotating-wave 
approximation is described by the Hamiltonian: (h = I; the waveguide group velocity is set to 
v,= l) 

H=J dkkalak + ViJ dyt [ajcri^^ + tr^^'W] + ^Qjcrf ^ 

+ ^Q2(rf + 72^ dk[alcr^^^ + (r^^^ak] (1) 
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where and Q2 are the atoms' transition frequencies; Vi^ are their respective coupling strengths 
to the waveguide field, and are related to the spontaneous decay rate of each atom by ^ = nV^ 
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Uf, (al) destroys (creates) a waveguide photon with energy k, and satisfies [ak,al,] = 6{k - k'). 

(12) 

a% ' are the lowering and raising operators for each atom. Once the S-matrix of the chiral mode 
is determined, the transport properties of the system in Fig.l may be obtained using standard 
techniques [3]. 

To solve the few-photon Fock-state transport properties for the Hamiltonian in Eq. (1), fol- 
lowing Ref. [24] we define the input, ain(0 = (2n)~^^^ J dke~^''^'~'°^ak(t[)), and output, aout(0 = 
(2n)~^^^ j dfce"'*^^'"''^a<:(?i), operators. Here, to and ti refer to times long before (to -00) and long 
after (ti ^ -1-00) the photons interact with the atoms. Following the procedure in [23, 24], we arrive 
at the input-output fomalism equations 

"W-i^^'-^^O, (2) 
r^^it) + i J^cr^\t)aM + "'"^2^'^'\ (3) 

(4) 

We note that the inclusion of a dipole-dipole interaction term in Eq. (l)oftheformg(cr*;.^V?^ -I- cr*^.^VL'^j, 
where g is the coupling rate, may be accounted for by making the replacement (tiTz)'^^^ ^ 
i'^iTiT^^^ + in the last term in Eqs. (3) and (4). The solution in this case would carry on in a 
similar fashion. Here, we exclude the dipole-dipole term. 

Below, we will solve Eqs. (3) and (4) to obtain the single -photon scattering amphtudes {p'\k^) = 
<0|aout(p)«J„(^)|0), and the two-photon scattering amplitude {pip~\kikp = {0\aout(Pi)aout(.P2)al^{k2)al^(ki)\0}. 
Here the k's and p's are incident and outgoing free -photon energy, respectively. One-Photon Scat- 
tering MatrixOne-Photon Scattering Matrix 



ni. ONE-PHOTON SCATTERING MATRIX 



Fourier transforming of Eq. (2) leads to the single-photon scattering amplitude 

{p-\k^) = <0|ai„(p)|r ) - i y J^(0\cr'f'\p)m (5) 
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FIG. 2. (Color online) Plot of the normalized resonance fluorescence, \B{ki^2, Pi,2)/t\^, assuming t\ = T2 = 
T = l/y. Q.C, Q-d and A,> are defined in the main text. Let E = (E - 2Q.c)t and A,_o = A,;oT. (a) E = 3, 
Q.J = y. (b) £ - 3, Orf - O.Sy. (c) £ = 3, = 0. (d) £ - 0. 



f 1 2) (12) 

In order to calculate (0|cr_' (p)|^+), we use Eqs. (3) and (4) and solve for {Q\(fl'^\tW)- The 

(12) (12) 

solution results in (0|cr_' \p)\k*) = s\ ' 6(k - p), where 

' VT:(fc-Q,+i-L)(fc-Q, + i-L) + _L. 

are the excitation amplitudes of the two atoms. We plug Eq. (6) into Eq. (5) to obtain the one-photon 
scattering matrix {p'\k^) = tk5{k - p) where 



^^_ (^-a-i^)(fc-Q.-i^) + 4 
'"(^-Qi+i^)(^-Q, + i^) + 4' 



IV. TWO-PHOTON SCATTERING MATRIX 

The two-photon scattering matrix can be written as 

(0|aout(iCl)<3out(iC2)«j'n(^2)fi'Jn(^l)|0) = </?tIam(P2) " i ^^^CT- ^Pl) " i ^^^CrL^\j!?2)l^l^2 > 

where we have used Eq. (5). We are then tasked with calculating the matrix elements (p||crl^'^^(p2)l^i^2 
For this purpose we again use Eqs. (3)-(4) to get a coupled differential equation for {p\\a'"l'^\t)\k\k1^) 
with inhomogenous terms of the type {p\\a^l'^\t)aiJj:)\kik2) and {p\\cr^p{t)(^^-^(f)\k\k2). From the 



definition of a^it) one can straightforwardly show that [24] 



ik2t<. 



In solving for <p+|crl'^(0cri^V)l^i^2 )' "^^ ^^^^^ 



+ 1^(2)/- 



As a result, we are left with calculating the matrix element {0\cr''_}\t)cr''^\t)\k[k2). From Eqs. (3) 
and (4), we can derive an operator equation 

^ [cr'y(t)cr'^(t)] = -i|Qi + - i^ - i^ j cri^^O^r^V) 



+ i J -cr'^^\t)(T^^\t)aM + i -x/ -crf\t)cr'^y{t)aUt), 



(7) 



where we have used the operator identity 



[aUt),(r'^!'^\t)]=0. 



(8) 



[Equation (8) is proved in Appendix A.] Equation (7) can then be used to solve for {0\cr^_^\t)(T^P(t)\kik2). 
At this stage, we have calculated all the prerequisite inhomogeneous terms in the coupled differential 
equations for (p||cr_' '(t)\kik2). The resulting two-photon S-matrix is 

{0KutiPi)aout(P2)al^(k2)al^(ki )|0) 
= tp.tp^ [6(ki - pi)6{k2 - pi) + 5{ki - P2)5{k2 - p\)\ 



+ { 



+ 



Pi Pi Pi Pi 

ny/TiT2 



(£,-Qi-Q2 + i^+i^) 
X 6{Ei - Eo) 



(9) 



where Ei = ki + k2 and Eo = pi + P2- 
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Knowing the S-matrix of the chiral model which contains the photon-atom interaction, the full 
S-matrix for the waveguide system in Fig. 1, which has both the left and right going photons, can 
then be constructed straightforwardly [3]. In particular, the fluorescence spectra of the transmitted 
and reflected photons are described by the same last three terms in Eq. (9). Below, we will discuss 
the results for the waveguide system shown in Fig.l. 

Examining Eq. (9), the first term describes an uncorrelated transport process where the en- 
ergy of individual photons is conserved. The following three terms, which we collectively label 
B{k\ 2, pi2)5{Ei - Eg), describe the fluorescence process where only the total energy of the photons, 
but not the individual energies, is conserved. In particular, the second and third terms represent 
fluorescence from each individual atom. The fourth term arises from the joint fluorescence in which 
both atoms are excited simultaneously, and contains a two-photon pole. These different fluorescent 
pathways interfere coherently, leading to a complex set of interesting efi'ects. 



V. RESULTS 

A. Fluorescence linewidth narrowing 

In Fig. 2 we plot the spectrum of fluorescence \B{ki^2, P\,2)/t\^, in the two-dimensional space 
spanned by the A, = (^i - k2)/2 and = (pi - P2)/2 axes. We assume the two atoms have 
identical waveguide coupling rates of y = l/r, but generally diflrerent resonant frequencies, and 
define Dc = "'2"^ ; i^d = "'2"^ - In Figs. 2(a)-2(c), we plot the fluorescent spectrum for incident 
two photons with total energy = 2Q.c + Sy. Examining Eq. (9), we see that in the A, - A^ 
plane, the poles of B{ki2, PU2) are the same as the poles of the atomic excitations (4i'f ' ■^pi?) with 
single photon input. The atomic excitation exhibits a sub-radiant state with poles corresponding 
to outgoing momenta pi 2 ~ O^c - i^^/T' and a super-radiant state with poles corresponding to 
outgoing momenta pi 2 = - 2iy. Consequently when the atoms' detuning is smaller than their 
linewidth (Q^ < y), the fluorescent features are quite narrow [Fig. 2(b)]. At zero detuning, however, 
the sub-radiant state has zero linewidth and no longer couples to externally incident photons. As 
a result, the fluorescence is dominated by the super-radiant poles [Fig. 2(c)], and the fluorescent 
hnewidth is doubled compared to a single atom's fluorescence. 
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FIG. 3. (Color online) P^^^{x), the joint detection probability density of the transmitted two-photon state. 

(a) k\ = Q.\, Q.ii - 2y, (^2 - ^^2) = {2, 1, in solid black, dotted blue and dashed red curves, respectively. 

(b) ki = Q.\, Q.d = 6y, (k2 - O.2) = {6, 0.25)y in solid black and dashed red curves, respectively. In 
(a) and (b), insets show the frequency kj overlaid on the single-photon transmission \ik\^. (c) ki = Q.i, 
k2 - O.C - {0.75, 0.5, 0.25)7 for the black (solid), blue (dotted) and red (dashed) curves, respectively. 



B. Fluorescence quenching and two-photon resonance 



Fluorescence is typically an unavoidable signature of interaction between an atom and multiple 
photons, since it arises from the inelastic scattering of one photon off an excited atom [3]. Here, as 
shown in Fig. 2(d), the fluorescence completely vanishes provided the total energy of the incident 
photons satisfies Et = Qi +Q.2- Fluorescence quenching was previously noted in a driven three-level 
system [31]. Here we show that quenching can occur in a system with two two-level systems. 

The effect of fluorescence quenching is closely related to the existence of a two-photon pole in the 
S-matrix, which provides the necessary pathway to cancel the contribution from the fluorescence 
of individual atoms. The existence of a two-photon pole indicates that the two atoms can be 
simultaneously excited by two photons, as long as the sum of the photon energy is near the sum of 
the transition energy of the two atoms. 

Whether a two-atom system can exhibit a two-photon pole has been an interesting question. In 
Ref. [32] it is argued that in the absence of dipole-dipole interaction [33], simultaneous excitation 
is not possible classically. The authors have also shown that simultaneous excitation is possible 
when the photon pair is frequency entangled, or when the two atoms interact via a quantized cavity 
field [34]. Our contribution is in showing the connection between the two-photon pole in the joint 
excitation of the atoms and fluorescence quenching. 
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C. Generation and control of two-photon bound states 



For an incoming state comprising two right-going photons with individual energies ki and k2, 
the resultant transmitted two-photon state is: 

Wr) = J dXidX2 [tkJkiSkiMi^l^^l) + lH(XuX2)] \Xi,X2)rr. 

Here 4 = (?/t-i-l)/2 is the single-photon transmission amplitude. S ki,k2(xi^X2) = ^e'(i'i^i+i'2X2) ^ Qi{kiX2+k2Xi)] 
and \xi,X2)rr = ^4(-^i)4(-^2)|0) where [cr(x),4(x')] = 6ix - x'). 



H(Xi,X2) = 



(1) 



2 4/ + h2 



n 



Fx{Xi,X2)\ \ - 



+ F2(XuX2) 



+ 



,(2) 



(2) 



2 % + ^2 

n 



Fi{Xi,X2) 



+ F2{Xi,X2)\l 



-1 
0^2 

i 



where 



V2e'^'( 



-X 



^l,2(-^l,^2) = „ „ 

Ti2DaDh 



Da+Di, 



D„-Du 



with D, 2 = (Da + Dh) 12, Da = Ei - 2a + i/Ti + i/T2, Dt, = [4Q2 + 41^/(1 /ti - I/T2) - (1/ti + 
1/^2)^] contains an uncorrelated-transport extended plane- wave term [ tkJk2S h^i^i^ ^2) \, 
and a bound-state term, X2)], which arises directly from fluorescence. 

Since X2|i/'«)p = {^R\a\{x2)a\{X])aR{x2)aR{xi)\il/R) is proportional to the joint-detection 

probability density Pj^^Xi - X2), this wavefunction therefore can be experimentally probed in a 
Hanbury Brown-Twiss coincidence measurement, with x = xi - X2 being the difference in optical 
path length from each detector to the beam splitter. 

The transmission amplitude 4 vanishes when k = 2, making it possible to eliminate the 
uncorrelated part of Ii/tr) by choosing ki = Q.i, which we do in Fig. 3. By doing so, the transmitted 
two-photon wavefunction is entirely described by the function H{xi,X2), which represents a two- 
photon bound state that decays with respect to the photon spacing x = xi - X2- 

When ^2 = ^c, P'f Kx) exhibits a bunching behavior with a global maximum at x = 0, as 
shown in Fig. 3(a). As one decreases ^2 from this value, a local minimum at .x; = starts to 
develop, indicating anti-bunching [35], which increases as k2 ^^2- this system, therefore. 



the two-photon bound-state can exhibit either a bunching or an anti-bunching behavior. This is 
in contrast to the one-atom case [3], where the bound state by itself is always bunched. It is also 
different from the typical resonance fluorescence experiments with classical input state, where 
anti-bunching is observed [36]. 

In Fig. 3(b) the transmitted bound state is plotted for a large atomic detuning of = 67, 
exhibiting spatial oscillation or quantum beats with a period of 2n/\Dh\. A two-photon plane wave 
S ki,k2(^i^^2) has a Pf'ixi - xj) that oscillates with a spatial period of k\ - ^2 where ki^ are the 
energy of the two individual photons. Similarly, the quantum beats here indicate that the bound state 
is two colored, with the single-photon energies centered approximately at Qi and 0.2- Previously, 
the second-order correlation function of two atoms driven by a coherent laser field was shown to 
display similar oscillations which are governed by the atomic detuning [37]. Here we note the 
connection between these oscillations and the spatial wavefunction of the two-photon bound state. 

Finally, we note that the spatial extent of the bound state is strongly dependent upon the atomic 
detuning Q^/, in the region where Q.^ < y. In Fig. 3(c), P^'jKx) is plotted for three values of Q.4 
satisfying O-dly < 1, with k2 - O.2 = O.^. For these values of Q.^, the sub-radiant poles of the 
scattering matrix, whose imaginary parts strongly depend on Q^, dominate the bound-state response. 
Consequently, as shown in Fig. 3(c), the spatial extent of the bound state can be significantly wider 
than the one-atom bound state. Taken together. Figs. 3(a)-(c) demonstrate the significant ability of 
the two-atom system to control the properties of the two-photon bound state, including its spatial 
extent, quantum beats and statistics. 

D. Experimental considerations 

We end by discussing some practical considerations relevant in experimental study of this system. 
All predicted effects in the paper require that the atomic resonance frequencies be close to each 
other. With respect to the tuning of atom resonance frequencies, we note that tuning of quantum dots 
in the optical frequency range has been achieved via magnetic fields [20], the application of a dc 
voltage [21], and through thermal heating [22]. In the microwave, the ability to independently tune 
the transition frequencies of individual qubits has recently been demonstrated [5]. It is conceivable 
that some of these techniques can be further developed to achieve independent tuning of two closely 
spaced qubits, particularly in the microwave frequency range. 

We also assumed that each atom predominantly couples to the waveguide, i.e. that the system 
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has a high jS factor. (The yS factor measures the fraction of the spontaneous emission going into the 
guided mode.) A high [5 factor has been reported in the experiments of Ref. [38, 39]. A value of 
jS < 1 may be accounted for by replacing each resonant frequency Q.12 with Q.12 - iyf2 where yj"!^ 
are the respective coupling rates of each atom into non-guided modes. Consequently, a non-unity 
jS will lead to the further narrowing of the fluorescence features in Fig. 2(b), and to the further 
broadening of the fluorescent features in Fig. 2 (a),(c). 

The main result of this paper [i.e., the two-photon S-matrix in Eq. (9)], is valid for two atoms 
with either identical or non-identical atom-waveguide coupling rates. As for the predicted efl'ects, 
linewidth narrowing and the various properties of the bound state persist for atoms with non- 
identical waveguide coupling rates. Complete fluorescence quenching, however, requires that the 
coupling rates be identical. For this purpose, we note that, in the microwave frequency range, Ref. 
[5] has demonstrated the capability of tuning the qubit coupling rate. 

We have calculated the response of the system to a two-photon Fock state input. The predic- 
tions about the properties of the transmitted two-photon state can be observed by a correlation 
measurement with a weak coherent state input, which can be generated with an attenuated laser 
beam [3]. The frequency linewidth of such a beam can in principle be made narrower than any 
of the spectral features that we predict here. Alternatively, we note the recent development of 
deterministic single-photon sources [40] as well as the demonstration of single photon pulses with 
arbitrary temporal shapes [41], both of which may facilitate the experimental study of this system. 

In our calculations, we have presented the S-matrix for two incident photons, each having 
well-defined energy ki and ^2- Any experiment, of course, uses a source with non-zero spectral 
bandwidth. The spectrum of the output state is a product of the scattering matrix and the spectrum 
of the input state, and thus can be directly calculated using Eq. (9). Since the fluorescence is 
already in the spectral domain, the use of a non-zero bandwidth input should not affect the results 
of such spectral measurement. In Fig. 3, we have shown that a pure two-photon bound state can be 
generated when one of the incident photons is on resonance with one of the atoms. The use of a 
photon pulse would therefore result in a background amplitude due to uncorrelated transport, in 
addition to the bound state. This background, however, can be made to be very weak, provided that 
the incident spectrum is significantly narrower than the atomic linewidth. 
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Appendix A: PROOF OF EQ. (8) 

From the Hamiltonian in Eq. (1), an equation of motion for ak(t) may be derived 

^at(t) = -iaM - iV:(r'^}\t) - iV2cr'-^\t) (Al) 
at 

Equation (Al) may be solved by integrating from time to: 

ak(t) = ak(to)e-"'^'''°^ - iVy f df'e'^^''" V^'\f') - iVj f dt'e"^''-"'W^^\t') (A2) 

Jto Jto 

Multiplying Eq. (A2) by and integrating over k while taking the limit to -oo, we obtain: 



aUt) = -^Jdk a,(t) + i + i 



where we have used the definition of ain(0 from Sec. II. It then follows that 



[aUt), cr^^\t)] = ^fdk [a,(t), cr^^^t)] + i ^ [cr^!\t), cr^^'^t)] + i ^ [cr^^(t), cr^}\t)] = 
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